Noise reduction method

ABSTRACT

A noise reduction method partitions frequency band into multiple sub-bands and estimates the signal-to-noise ratio (SNR) value for each sub-band. An over-subtraction factor of each sub-band is determined based on the estimated SNR value. Then. the clean speech spectrum estimate is determined by performing spectral over-subtraction on each sub-band, so as to determine the clean speech signal from the estimated clean speech spectrum.

BACKGROUND OF THE INVENTION

[0001] 1. Field of the Invention

[0002] The present invention relates to a noise reduction method and, more particularly, to a method using spectral subtraction to reduce noise.

[0003] 2. Description of Related Art

[0004] The spectral subtraction method has been proven effective in enhancing speech degraded by additive noise. It is simple to implement, hence is suitable as the pre-processing scheme for speech coding and recognition applications. This method subtracts the noise spectrum estimate from the noisy speech spectrum to estimate the speech magnitude spectrum, so as to obtain the clean speech signals.

[0005]FIG. 1 shows the flowchart of the aforementioned spectral subtraction method, wherein the input noisy speech is divided into a plurality of continuous frames, and each frame is represented by an additive noise model:

y _(r)(k)=s _(r)(k)+w _(r)(k),

[0006] where y_(r)(k), s_(r)(k) and w_(r)(k) denote respectively the k-th noisy speech, clean speech, and noise sample of the r-th frame. Taking the fast Fourier transform of the noisy speech frame y_(r)(k) (step S101), the noisy speech spectrum of the r-th frame at the k-th frequency component is obtained and denoted as |Y_(r)(k)|². In addition, the noisy speech y_(r)(k) is also applied in a silence detection process (step S102) and a noise spectrum estimation process (step S103) to estimate a noise spectrum, denoted as |W_(r)(k)|². After performing a spectral subtraction process (step S104), the energy spectrum of clean speech is obtained as follows:

|Ŝ _(r)(k)|² =|Y _(r)(k)|² −|W _(r)(k)|²  (1)

[0007] If the phase spectrum of the clean speech can be approximated by the phase spectrum of the noisy speech, the estimate of clean speech ŝ_(r)(k) can be obtained by taking the inverse fast Fourier transform of |Ŝ_(r)(k)|².

[0008] Such a method is suitable as the pre-processing scheme for speech coding and recognition applications because it is easy, effective and simple to implement. However, the noise spectrum estimate may cause a relatively large spectral excursion in the spectrum estimate of clean speech. This spectral excursion will be perceived as time varying tones contributing to the so-called musical noise.

[0009] To reduce the musical noise Berouti et al proposed a noise reduction method to over-subtract the noise spectrum estimate, and a description of such can be found in M. Berouti, R. Schwartz, and J. Makhoul “Enhancement of speech corrupted by acoustic noise”, pp. 208-211, 1979 IEEE, which is incorporated herein for reference, wherein the formula (1) is modified as:

|Ŝ _(r)(k)|² =|Y _(r)(k)|²−α_(r) ·|W _(r)(k)|²·α_(r)≧1  (2)

[0010] so as to decrease the influence caused by the excursion of the noise spectrum estimate and thus reduce the effect of musical noise. In the method, the over-subtraction factor α_(r) was determined by the signal-to-noise ratio (SNR) of the processing frame, and can be expressed by formula: $\begin{matrix} {{\alpha_{r} = {\alpha_{0} + {{SNR}_{r} \cdot \frac{1 - \alpha_{0}}{{SNR}_{1}}}}},} & (3) \end{matrix}$

[0011] where α₀ is pre-selected over-subtraction factor when SNR=0, SNR₁ is pre-selected SNR value when α_(r)=1, SNR_(r) is the estimate of signal-to-noise ratio of the processed r-th frame. Based on the formula (3), it is known that α_(r) is inversely proportional to SNR_(r). The smaller the SNR_(r) is, the larger the α_(r) is, and a larger α_(r) is helpful in removing the larger noise spectrum excursion.

[0012] Examining human speech spectrum, it is known that the speech energy distributes non-uniformly and often concentrates on lower frequency components. Hence SNR differs with frequencies and often have larger values at lower frequency components. From the formula (3), it is known that more suppression is needed for lower SNR and vise versa. High-frequency components thus need more suppression to avoid musical noise, while low-frequency components need less suppression to prevent speech distortion. However, for the over-subtraction method based on formulas (2) and (3), it faces the problem of too much over-subtraction and hence speech distortion at low-frequency components while too less over-subtraction and hence musical noise at high-frequency components. Accordingly, improved schemes are proposed to avoid such a problem, and one of the schemes can be found in Kuo-Guan Wu and Po-Cheng Chen “Efficient speech enhancement using spectral subtraction for car hands-free application”. 2001 Digest of technical papers, pp. 220-221, which is incorporated herein for reference. However, it is unable to completely eliminate the problem. Therefore, there is a need for the above conventional noise reduction method to be improved.

SUMMARY OF THE INVENTION

[0013] The object of the present invention is to provide a noise reduction method capable of effectively eliminating the musical noise and reducing speech distortion.

[0014] To achieve the object, the noise reduction method divides input noise speech into a plurality of continuous frames, determines noisy speech spectrum for each frame, and partitions frequency band into multiple sub-bands to determine clean speech spectrum from the noisy speech spectrum on each sub-band. The method is provided to first estimate noise spectrum of r-th frame at k-th frequency component from the noisy speech of r-th frame by silence detection and noise spectrum estimation. Next, the signal-to-noise ratio (SNR) value of i-th sub-band for r-th frame is estimated. Then, an over-subtraction factor of sub-band i is determined based on the estimated sub-band SNR. Finally, the clean speech spectrum estimate is determined by performing a spectral subtraction on each sub-band.

[0015] Other objects, advantages, and novel features of the invention will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings.

BRIEF DESCRIPTION OF THE DRAWINGS

[0016]FIG. 1 is the flowchart of a conventional spectral subtraction method.

[0017]FIG. 2 is the flowchart of the noise reduction method in accordance with the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

[0018] With reference to FIG. 2, there is shown the flowchart of a preferred embodiment of the noise reduction method in accordance with the present invention. As shown, the input noisy speech of the r-th frame y_(r)(k)=s_(r)(k)+w_(r)(k) is processes by FFT (fast Fourier Transform) (step S201) to obtain its energy spectrum |Y_(r)(k)|². The noisy speech y_(r)(k) is also processed by silence detection (step S202) and noise spectrum estimation (step S203) to estimate the noise spectrum of the r-th frame, denoted as |W_(r)(k)|².

[0019] For the noisy speech spectrum |Y_(r)(k)|² and noise spectrum |W_(r)(k)|², the method of the present invention utilizes a sub-band over-subtraction mechanism to determine the estimate of clean speech spectrum |Ŝ_(r)(k)|², which is then processed by IFFT (Inverse Fast Fourier Transform) (step S207) for being restored to enhanced frame signal ŝ_(r)(k). The method of the present invention partitions the frequency band into multiple sub-bands and perform over-subtraction on each sub-band. To implement over-subtraction on each sub-band, it is first performed a sub-band SNR estimation (step S204) to estimate a SNR value for determining the over-subtraction factor of the sub-band. The SNR value can be obtained by a regression formula as follows: ${{{SNR}_{r}(i)} = {{\mu \cdot {{SNR}_{r - 1}^{o}(i)}} + {\left( {1 - \mu} \right) \cdot 10 \cdot {\log_{10}\left( {\frac{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{Y_{r}\left( {i,k} \right)}}^{2}}{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{W_{r}\left( {i,k} \right)}}^{2}} - 1} \right)}}}},$

[0020] where i is the index of sub-band, SNR_(r)(i) is SNR estimate of the i-th sub-band for the r-th frame, |Y_(r)(i,k)|² is noisy speech spectrum of the r-th frame at the k-the frequency component of the i-th sub-band, |W_(r)(i,k)|² is the corresponding noise spectrum, 0<μ<1, and SNR_(r−1) ^(o)(i) is the SNR of the sub-band for the previous frame after noise reduction, which is expressed by the following formula: ${{{SNR}_{r - 1}^{o}(i)} = {{10 \cdot \log_{10}}\frac{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{{\hat{S}}_{r - 1}\left( {i,k} \right)}}^{2}}{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{W_{r - 1}\left( {i,k} \right)}}^{2}}}},$

[0021] where |Ŝr−1(k,i)|² is the estimate of clean speech spectrum of the previous, i.e., the (r−1)-th, frame after being processed in the sub-band i.

[0022] In step S205, the sub-band over-subtraction factor α_(r)(i) is determined based on the estimated sub-band SNR value SNR_(r)(i), and is expressed by the formula as follows: ${{\alpha_{r}(i)} = {{\alpha_{0}(i)} + {{{SNR}_{i}(i)} \cdot \frac{1 - {\alpha_{0}(i)}}{{SNR}_{1}(i)}}}},$

[0023] where α₀(i) is pre-selected over-subtraction factor when the actual SNR_(r)(i)=0 at sub-band i, and SNR₁(i) represents pre-selected SNR value when α_(r)(i)=1.

[0024] Once determining the over-subtraction factor α_(r)(i) for each sub-band i, it is able to perform spectral over-subtraction on each sub-band i (step S206), as expressed by the following formula:

|Ŝ _(r)(i,k)|² =|Y _(r)(i,k)|²−α_(r)(i)·|W _(r)(i,k)|²,

[0025] wherein the determined |Ŝ_(r)(i,k)|² is the clean speech spectrum at sub-band i for the r-th frame. After performing over-subtraction for each sub-band i, the IFFT is applied (step S207) to obtain the estimated enhanced frame signal ŝ_(r)(k).

[0026] In executing the aforementioned method, due to the small number of frequency samples in the lower bands, there will be large variation in sub-band SNR estimate when the noise is strong, which may cause an error in α_(r)(i) and influence the quality of the restored speech. To avoid such a problem, in step S205, the SNR value SNR_(r) of the whole frame is incorporated into modification of sub-band over-subtraction factors as follows:

[0027] α_(r)(i)=α_(max) if SNR_(r)<SNR_(min),

[0028] where SNR_(min) is pre-selected minimum value of SNR.

[0029] Furthermore, in this embodiment, the step S204 employs regression scheme to estimate the SNR value for determining the over-subtraction factor of the sub-band. However, in practical application, the SNR value of sub-band can also be determined by other known speech signal SNR estimation methods, for example, the high order statistic method described in Elias Nemer, Rafik Goubran and Samy Mahmoud: ‘SNR estimation of speech signals using subbands and fourth-order statistics’, IEEE Signal Processing Letters, 1999, vol. 6, no. 7, pp. 171-174, which is incorporated herein for reference.

[0030] To verify the effect of the present noise reduction method, noisy speech data is generated by adding clean speech data with white Gaussian noise of variant magnitudes to form 3 segmental SNRs: 15 dB, 10 dB and 5 dB. Eight clean speech sentences are collected with 5 sentences from males and 3 from females. Table 1 compares the averaged segmental SNR improvements of conventional over-subtraction method (with parameters of α₀=7.5 and SNR₁=20) and those of the present method (with parameters of α₀(1˜18)=2, SNR₁(1˜13)=1.5, SNR₁(14˜18)=1.25) with sub-band SNR obtained from clean speech data. TABLE 1 Method Present Conventional sub-band Improvement of Input SNR over-subtraction over-subtraction the present method 15 dB 2.39 3.33 39.3% 10 dB 3.86 4.76 23.3%  5 dB 5.64 6.64 17.5%

[0031] From this comparison, it is known that at 15 dB input SNR, the present method has the potential of achieving 40% improvement over the conventional method. The potential improvements increase with input SNR.

[0032] Table 2 compares the averaged segmental SNR improvements of conventional over-subtraction method (with parameters of α₀=7.5 and SNR₁=20) and those of the present method (with parameters of α₀(1˜18)=2, μ=0.25, SNR₁(1˜9)=10, SNR₁(10˜13)=15, SNR₁(14˜16)=2, and SNR₁(17˜18)=1.25) with sub-band SNR obtained from the step S204 of sub-band SNR estimation. TABLE 2 Method Present Conventional sub-band Improvement of Input SNR over-subtraction over-subtraction the present method 15 dB 2.39 2.80 17.0% 10 dB 3.86 4.09  6.0%  5 dB 5.64 5.96  5.7%

[0033] From Table 2, it is known that at input SNR=15 dB, although the SNR value of sub-band is obtained by estimation, the present method still can achieve 17% improvement over the conventional method.

[0034] Although the present invention has been explained in relation to its preferred embodiment, it is to be understood that many other possible modifications and variations can be made without departing from the spirit and scope of the invention as hereinafter claimed. 

What is claimed is:
 1. A noise reduction method for dividing input noise speech into a plurality of continuous frames, determining noisy speech spectrum for each frame, and partitioning frequency band into multiple sub-bands to determine clean speech spectrum from the noisy speech spectrum on each sub-band, the method comprising: (A) estimating noise spectrum |W_(r)(k)|² of r-th frame at k-th frequency component from the noisy speech y_(r)(k) of r-th frame by silence detection and noise spectrum estimation; (B) estimating signal-to-noise ratio (SNR) value SNR_(r)(i) of i-th sub-band for r-th frame; (C) determining an over-subtraction factor α_(r)(i) of sub-band i based on the estimated SNR_(r)(i); and (D) determining clean speech spectrum estimate by performing, on each sub-band, a spectral subtraction: |Ŝ _(r)(i,k)|² =|Y _(r)(i,k)|²−α_(r)(i)·|W _(r)(i,k)|², where |Y_(r)(i,k)|² is noisy speech spectrum of the r-th frame at the k-the frequency component of the i-th sub-band, |W_(r)(i,k)|² is corresponding noise spectrum, and |Ŝ_(r)(i,k)|² is clean speech spectrum at sub-band i for the r-th frame.
 2. The noise reduction method as claimed in claim 1, wherein in step (C), the over-subtraction factor of the i-th sub-band for the r-th frame is: ${{\alpha_{r}(i)} = {{\alpha_{0}(i)} + {{{SNR}_{r}(i)} \cdot \frac{1 - {\alpha_{0}(i)}}{{SNR}_{1}(i)}}}},$

where α₀(i) is pre-selected over-subtraction factor when the actual SNR_(r)(i)=0 at sub-band i, SNR₁(i) represents pre-selected SNR value when α_(r)(i)=1, and SNR_(r)(i) is SNR estimate of the i-th sub-band for the r-th frame.
 3. The noise reduction method as claimed in claim 2, wherein, the over-subtraction factor α_(r)(i) of the sub-band is modified by the SNR value SNR_(r) of the frame as: α_(r)(i)=α_(max) if SNR_(r)<SNR_(min), where SNR_(min) is a pre-selected minimum value of SNR.
 4. The noise reduction method as claimed in claim 2, wherein the SNR_(r)(i) is obtained by a regression process: ${{{SNR}_{r}(i)} = {{\mu \cdot {{SNR}_{r - 1}^{o}(i)}} + {\left( {1 - \mu} \right) \cdot 10 \cdot {\log_{10}\left( {\frac{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{Y_{r}\left( {i,k} \right)}}^{2}}{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{W_{r}\left( {i,k} \right)}}^{2}} - 1} \right)}}}},$

where 0<μ<1, and SNR_(r−1) ^(o)(i) is the SNR of the sub-band i for the previous frame after noise reduction.
 5. The noise reduction method as claimed in claim 4, wherein SNR_(r−1) ^(o)(i) is determined by: ${{SNR}_{r - 1}^{o}(i)} = {{10 \cdot \log_{10}}{\frac{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{{\hat{S}}_{r - 1}\left( {i,k} \right)}}^{2}}{\sum\limits_{k \in {{sub} - {{band}\quad i}}}{{W_{r - 1}\left( {i,k} \right)}}^{2}}.}}$


6. The noise reduction method as claimed in claim 2, wherein the SNR_(r)(i) is obtained by a high order statistic method.
 7. The noise reduction method as claimed in claim 1, wherein the noisy speech is processed by fast Fourier transform to obtain the noisy speech spectrum.
 8. The noise reduction method as claimed in claim 1, wherein the noisy speech is processed by silence detection and noise spectrum estimation to estimate the noise spectrum.
 9. The noise reduction method as claimed in claim 1, wherein in step (D), the determined clean speech spectrum estimate is processed by inverse fast Fourier transform to obtain corresponding enhanced speech signal. 